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SOME  BASIC  CONCEPTS  FOR  MAGNET  COIL  DESIGN 

W.  F,  Gauster* 


In  the  United  States  and  abroad  there  is  intensive  activity  in  thermonuclear  research J 
A  major  part  of  this  effort  is  directed  toward  containment,  in  strong  magnetic  fields,  of 
plasma  of  the  necessary  high  particle  energies,  corresponding  to  temperatures  of  hundreds 
of  millions  of  degrees  Kelvin.  Magnetic  field  design  in  thermonuclear  research  has  already 
been  developed  to  a  relatively  broad  field  of  electrical  engineering.^  The  purpose  of  this 
paper  is  to  discuss  some  theorems  which  are  of  special  interest  for  the  design  of  large  d-c 
magnet  coils  as  used,  for  instance,  in  the  DCX  (Direct  Current  Experiment)  thermonuclear 
research  program  of  the  Oak  Ridge  National  Laboratory.^ 

I.  CYLINDRICAL  MAGNET  COILS  WITH  SQUARE  ENDS  AND 
UNIFORM  CURRENT  DENSITY 

The  design  of  high  power  magnet  coils  Is  based  on  fundamental  work  by  Fabry, 
Cockcroft,^  and  Bitter.^  It  is  advisable  to  summarize  a  few  well-known  points  of  this 
theory.  The  simplest  type  of  magnet  coil  is  cylindrical  with  square  ends  and  is  designed 
for  operation  with  uniform  current  density  i  (Fig.  1).  Fabry  showed  that  the  magnetic 
field  strength  at  the  coil  center  can  be  expressed  by: 


where 

P  =  ohmic  loss  power,  in  watts, 

conducting  cross  section 

k  =  “space  factor^*  - ; - - - ; - ; — ; —  , 

total  cross  section  including  insulation 

*Oak  Ridge  National  Laboratory,  operated  by  Union  Carbide  Corporation  for  the  U.S.  Atomic 
Energy  Commission.  The  author  wishes  to  express  his  gratitude  to  Dr.  J.  P.  Neal  and  Mr.  T.  F, 
Connolly  tor  their  helpful  co-operation  in  editing  this  paper. 
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p  =  resistivity  of  the  coil  conductor,  in  ohm-cm, 

=  inside  radius  of  the  coil,  in  cm. 

//q  Is  measured  in  oersteds,  and  Bq  has  the  same  numerical  value  In  gauss  (no  ferromag¬ 
netic  material  is  supposed  to  be  in  the  magnetic  field).  G  is  a  dimensionless  factor,  the 
so-called  "Fabry  factor,"  which  does  not  depend  on  the  size  but  only  on  the  shape  of  the 
coil. 

By  introducing  the  ratios: 


h 

a. 


(2) 


(see  Fig.  1),  G  can  be  expressed  as  a  function  of  a  and  in  an  elementary  way.^°  It  is 
convenient  to  represent  G  by  a  family  of  curves  (see  Fig.  2).^  A  "flat"  maximum  of 
G  =  0.179  occurs  for  3.09,  ^  =  1.88. 

Fabry  also  used  a  "volume  factor  v/*  defined  by: 


total  coil  volume 


V  - 


a 


3 

1 


(3) 


For  instance,  in  the  case  of  a  cylindrical  coil  with  square  ends,  v  -  27t(o?  -  l)/3.  In 
Fig.  2f  V  -  constant  curves  are  shown  by  dashed  lines.  G^^^  corresponds  \o  v  ~  101, 
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Fabry  and  Volume  Factor  Diagrams  of  Cylindrical  Magnet  Coils  with 
Square  Ends  and  Uniform  Current  Density, 


Figure  2. 
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A  somewhat  smaller  G  can  be  achieved  with  a  much  smaller  volume  factor.  For  instance, 
the  DCX  coils  are  designed  with  the  following  data:  ® 

a  =  1.53, 

i8  =  0.618, 

V  =  5.21, 

G  =  0.128. 

Fabry  gave  a  table  showing,  for  volume  factors  between  1  and  10^,  values  of  a  and 
|S  for  which  G  assumes  maximal  values.^  The  DCX  coil  is  such  a  minimum-volume  coil. 

It  should  be  pointed  out,  however,  that  for  many  applications  a  minimum  volume  factor  is 
not  of  special  importance,  since  the  whole  field  configuration  produced  by  the  magnet  coil 
must  be  considered. 


n.  FABRY  FACTOR  OF  A  LAYER  COIL 


A  coil  shall  be  called  a  "layer  coil**  if  it  is  wound  in  straight  homogeneous  layers; 
i.e.,  if  current  density  i  (measured  in  amp/cm^)^ resistivity  p,  and  space  factor  X  are 
functions  only  of  the  distance  from  the  axis  of  r.  If  these  quantities  are  functions  of  r 
and  2,  the  coil  will  be  called  a  “loop  coil.**  Figure  3  shows  a  layer  coil  which  is  sym¬ 
metrical  to  a  plane  through  0  perpendicular  to  the  2r-axis.  In  the  following,  a  general  ex¬ 
pression  for  the  Fabry  factor  G  will  be  derived  for  a  symmetric  layer  coil  where  the 
strength  is  referred  to  the  coil  center.  A  similar  expression  can  be  found  for  an  asym¬ 
metric  layer  coil,  and  Hq  referred  to  any  point  at  the  axis  (see  Appendix  II). 

A  cylindrical  current  sheet  with  a  radius  r,  the  length  2z,  and  a  thickness  dr  produces 
in  the  center  the  field  strength: 


An 


An 


dHf.  = —  Xi  cos  4>  dr  ^  — Xi  -—dr  . 
0  10  10  s 


Therefore, 


An 


lo  J 


2 
r^a 


Xi--dr  . 
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The  power  dissipated  in  the  infinitesimal  current  sheet  is: 


(4) 


Therefore, 


dP  =  pz^  dV ^  ^  =  pXi^  dV  =  AnpXi^rz  dr  . 


P  =  An  j*  ^  pXi^rz 


dr 


(5) 


Normalized  (dimensionless)  quantities,  which  are  in  general  functions  of  r,  are  introduced 
in  the  following  way: 


.+ 

i 


(6) 


The  subscript  1  refers  to  the  innermost  layer  where  r  =  a 
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(4a) 


Optimum  Current  Density  of  a 
Symmetrical  Layer  Coil 
(Shape,  and  Known). 


‘^OPT-7^  Fig.  3. 

Figures  3  and  4. 
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(p‘^=  1)  will  be  considered.  The  current  density  is  also  assumed  to  be  constant,  in¬ 
versely  proportional  to  the  radius  r,  or  inversely  proportional  to 

a .  z  ^  =  1  . 

yj  377  /~^  CC—  1 

5  Vi  +  .Jo?  -  1 


The  maximum  value  for  G  is  0.172,  for  an  angle  =  45°  {k  =  1),  and  a=  2.7.  This 
is  case  4  of  Bitter's  paper. ^ 


G  - 


5  y  1  +  ^2  y  a  - 1 


(9) 


The  maximum  value  for  G  is  0.201,  again  for  an  angle  0  =  45°,  but  for  a=  4.5.  This 
is  Bitter's  case  5. 


The  maximum  value  for  G  is  0.250,  for  an  angle  ^  =  45°,  and  a  approaching  infinity. 
It  seems  that  this  case  has  not  been  discussed  previously;  however,  as  will  be 
shown,  it  deserves  interest  because  this  current  distribution  yields  the  maximum 
G  for  a  cylindrical  layer  coil  with  tapered  ends. 


IV.  OPTIMUM  CURRENT  DENSITY  DISTRIBUTION  OF  A  LAYER  COIL 

In  the  previous  paragraph  maximum  values  of  G  were  considered  which  occur  for 
special  values  of  k  and  a,  i.e.,  for  special  shapes  of  cylindrical  coils  with  tapered  ends. 
Another  problem  is  to  find  for  a  coil  of  any  definite  shape  a  current  density  distribution 
z(r)  at  which  G  becomes  a  maximum  (optimum  current  density  distribution).  Special  cases 
of  this  problem  were  discussed  by  Maxwell,  Fabry,  and  Bitter.^  Here  a  general  law 

for  the  optimum  current  density  distribution  of  a  coil  of  known  shape  and  a  general  ex¬ 
pression  for  the  Fabry  factor  (optimum  Fabry  factor)  which  results  in  this  case  are  de¬ 
rived.  For  simplicity  a  symmetrical  layer  coil  is  supposed,  and  is  referred  to  the  coil 
center.  The  more  general  problem,  to  find  optimum  current  density  distribution  and  op¬ 
timum  Fabry  factor  for  a  layer  coil  of  any  shape,  can  be  solved  in  a  similar  way.  As 
shown  in  the  Appendix  I,  a  certain  value  of  (Eq.  4)  can  be  produced  with  a  minimum 
power  P  (Eq.  5)  if  the  following  relation  is  satisfied. 

prsz  =  p ^  ^  =  constant.  (11a) 


5 


Thus 


.+ 

i 


1  ^1^1 


(nb) 


That  is  to  say:  a  desired  magnetic  field  strength  in  the  center  of  a  symmetrical  layer 
coil  is  achieved  with  minimum  power  when  the  product  of  specific  resistance,  radial  dis^ 
tance  of  a  layer  from  the  axis,  distance  of  the  layer  end-point  from  the  coil  center,  and 
current  density  are  constant. 

Substituting  z'*'  from  Eq.  11b  into  Eq.  7b,  we  obtain  the  optimum  Fabry  factor: 


G 


opt 


A" 

+ 

9 


z 


dr  (layer  coil). 


(12) 


Some  examples  will  be  discussed  below. 


V.  CURRENT  DENSITY  FACTOR  J  OF  A  LAYER  COIL 

As  shown  above,  the  Fabry  factor  G  is  dimensionless  and  characteristic  of  a  certain 
coil  shape  and  the  distribution  of  relative  current  density,  space  factor,  and  resistivity. 
The  absolute  value  of  is  found  by  multiplying  G  by  a  scaling  factor: 


Pl^l 

(see  Eq.  7a).  A  similar  procedure  may  now  be  used  to  find  the  absolute  value  of  the 
current  density  iy  From  Eq.  5a, 


47rpiA,  f  p\'^{i'^)^rz  dr 


(13) 


This  equation  can  be  written  as: 


(14a) 


(14b) 


The  dimensionless  factor  /  shall  be  called  the  “current  density  factor.’* 
value  of  iy  is  found  by  multiplying  /  by  a  scaling  factor: 


Pi  ^1 ^1 


The  absolute 


The  analogy  to  the  Fabry  factor  is  obvious. 
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Two  special  cases  might  be  considered  first.  For  a  coil  of  any  shape  with  constant 
resistivity,  space  factor,  and  current  density,  the  current  density  factor  becomes  simply 
(see  Eq.  3): 

/=_L  (p"  =  A"=f"=  1)  .  (15) 

In  this  special  case  the  current  density  factor  is  the  reciprocal  of  the  square  root  of  the 
volume  factor. 

For  a  symmetrical  layer  coil  of  any  shape,  operated  with  optimum  current  distribution 
(Eq.  nb),  the  current  density  factor  (Eq.  14b)  becomes: 


Or,  considering  Eq.  12: 


1 

j  - -  -  (layer  coil)  . 

105,  G  ' 

1  opt 


When  the  Fabry  factor  of  a  symmetrical  coil  with  optimum  current  distribution  is  known, 
the  absolute  values  of  the  current  densities  can  be  found  easily  by  Eqs.  17,  14a,  and  11b. 
For  a  coil  section  similar  to  that  shown  in  Fig.  3,  =  5^,  and  therefore 


1  1 

^opt=^^ - (layer  coil),  =S,  . 

Opt 


VI.  EXAMPLES  OF  LAYER  COILS  WITH  OPTIMUM  CURRENT 
DENSITY  DISTRIBUTION 


a.  Cylindrical  Coil  with  Square  Ends  and  Optimum  Radial  Current  Distribution, 
p''=\‘'=  1  (see  Fig.  1).'“’-^ 

Equation  11b  becomes  in  this  case: 


and  Eq.  12  yields: 


1  IItt  I  CL^(1  +  ^^) 


(19) 
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For  a  coil  with  a~  3,  (3  =  2,  the  Fabry  factor  is  =  0.197.  From  Eq.  17  the  current 
density  factor  is: 


} 


opt 


1  1 


0.226  . 


Equation  lib  yields  for  the  current  density  at  r  = 


a,s, 


0.206 z,  . 


It  is  interesting  to  compare  these  data  with  those  for  a  cylindrical  coil  with  square 
ends  and  constant  current  density  and  with  the  same  Uj,  a,  /3,  p,  and  X.  The  field 
strength  at  the  coil  center  in  both  coses  is  the  same.  Then  the  power  ratio  is  (see  Eq.  1): 


P 

~P 


/0.-179Y 
\  0-197/ 


0.825  . 


That  is  to  say,  17.5%  of  the  power  can  be  saved  by  optimizing  the  current  density  distri¬ 
bution.  From  Eq.  15: 


1 

/'=  —=0.100  . 

Therefore  (see  Eq.  14a), 


JG'  0.226  x  0.179  ^ 

i\~  j'^'  }'G  0.100x0.197 

The  current  density  in  the  optimized  coil  a\  r  -  is  2.05  iy  and  decreases  to  2.05 z|  x 

0.206  =  0.424  /;  at  r  = 


b.  Cylindrical  Coil  with  Tapered  Ends  and  Optimum  Radial  Current  Distribution, 
p ^  ^  A.^  =  1  (see  Fig.  4). 

Equation  lib  becomes  in  this  case: 


.+ 

I  = 


(20) 


A  cylindrical  coil  with  tapered  ends  (resistivity  and  space  factor  constant)  has  a  maximum 
Fabry  factor  when  the  current  density  is  inversely  proportional  to  the  square  of  the  dis¬ 
tance  from  the  axis.  This  is  exactly  the  previously  considered  third  example  of  Section 
III.  For  k  =  1  and  3,  Eq.  10  yields  =  0.204.  From  Eq.  17,  } =  0.346. 

From  Eq.  8  it  follows  that  for  a  cylindrical  coil  with  tapered  ends,  with  k  =  1  and 
a=  3  but  with  constant  current  density,  G'=  0.170.  With  the  same  in  each  coil, 

P_  /0.170^2 
"  V  0.2"04 


0.695 


8 


that  is,  30.5%  power  can  be  saved  by  optimizing  the  current  density  distribution.  From 
Eq.  15, 


J'=  - =  ■  =  0-0958  . 

Therefore, 

0.346  x  0.170 

—  =  -  =  3.01  . 

z(  0.0458  x  0.204 

The  current  density  in  the  optimized  coil  oi  r  =  is  3.05  and  decreases  to  of  this 
value,  i.e.,  to  0.334  z  j  at  r  =  a2*  Whether  the  absolute  value  of  the  maximum  current 
density  is  acceptable  or  not  depends  on  the  numerical  value  of  the  scaling  factor, 


VII.  SHAPE  OPTIMIZATION  OF  LAYER  COILS 

Section  IV  dealt  with  finding  the  optimum  current  density  distribution  of  a  layer  coil. 
Another  problem  is  that  of  the  optimum  shape  of  a  coil  (i.e.,  where  G  is  a  maximum)  when 
the  current  density  distribution  z(r)  is  known.  For  simplicity,  only  symmetrical  layer 
coils  will  be  considered.  The  case  of  an  unsymmetrical  layer  coil  is  discussed  in  Ap¬ 
pendix  II. 

It  can  be  shown  (see  Appendix  I)  that  the  optimization  condition  in  this  case  is: 


z  p(r^  +  z 

r  r 

We  consider  the  following  special  cases: 


(21) 


a.  p'^^  1  (Constant  Current  Density) 

Using  the  coordinates  5  and  0  (see  Fig.  3),  we  obtain  for  the  contour  curve  of  the 
coil  the  equation: 


— ; - =  constant.  (22) 

sin 

This  case  was  considered  by  Fabry. He  found  G  =  0.18. 

*  '  max 

b,  =  1,  z^  =  a^/r  (Current  Density  Inversely  Proportional  to  the  Distance  from  the  Axis) 

The  equation  of  the  contour  curve  is: 

s 

-  =  constant.  (23) 

sin^  ^ 

The  Fabry  factor  is  again  near  0.18. 
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c.  =  ]^  (a^/r)^  (Current  Density  Inversely  Proportional  to  the  Square  of  the  Dis¬ 

tance  from  the  Axis) 

In  this  case: 


(h  =  constant  =  -- 
^  4 


(24) 


This  layer  coil  is  a  coil  with  tapered  ends,  as  considered  in  sections  lllc  and  Vl^.  It 
can  be  shown  that  this  typ6  of  a  layer  coil  represents  the  solution  of  the  problem  of 
finding  the  optimum  Fabry  factor  when  both  shape  and  current  density  distribution  are 
optimized  simultaneously.  However,  in  this  case  the  maximum  current  density  is  rela¬ 
tively  very  high,  and  this  fact,  combined  with  the  difficulties  of  building  such  coils, 
seems  to  limit  their  practical  application. 


VIII.  LOOP  COILS 

As  mentioned  previously,  a  loop  coil,  in  contrast  to  a  layer  coil,  is  one  in  which  the 
quantities  z,  r,  p,  and  A  are  functions  of  both  coordinates  r  and  z.  Loop  coils  have  little 
practical  importance  (although  a  slightly  higher  G  is  attainable)  since  their  design  is  dif¬ 
ficult.  However,  some  basic  theorems  concerning  loop  coils  and  some  numerical  examples 
may  be  of  interest. 

Assume  that  a  loop  has  an  infinitesimally  small  cross  sectional  area  dA,  in  the  z-r 
plane,  and  that  its  location  and  radius  is  determined  by  the  coordinates  a  and  r  (distance 
from  the  origin  5  =  +  r^);  then  the  magnetic  field  strength  produced  at  the  origin  is:^® 

H  f K  —  i  JA  .  (25) 

“  10  •’a  ,3 

When  normalized  with  respect  to  an  appropriately  selected  reference  loop  (indicated  by 
the  subscript  1),  this  equation  becomes: 


277 

=  — X 
°  10 


dA 


A  s- 


The  power  dissipated  is: 


or,  normalized. 


P  =  277  C  pXrP-  dA  , 

^A 


P  =  27rp,X,f2  ^  _ 


Correspon(Jing  to  Eqs.  7a  and  7b  we  obtain  now  for  a  loop  coil: 


fPX, 


Wo  =  G 


p.a 


(25a) 


(26) 


(26a) 


(27a) 
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■\J  2ff<ai 


f  \^—i*dA 

rj- 

I  p*k^r{i*)^  dA 


(loop  coil) 


For  the  current  density  factor  we  obtain  (compare  Eqs.  14a  and  14b): 

/~P 


r(z'^)2  dA 


(loop  coil)  . 


J.  C.  Maxwell  solved  the  problem  of  designing  a  galvanometer  in  which  a  given  small 
electromotive  force  would  produce  the  greatest  possible  deflection.  ^  ^  In  a  more  general 
way  we  can  ask  about  the  optimum  current  density  distribution  of  a  loop  coil  of  any  shape, 
for  which  p  and  X  are  known  functions  of  ^  and  r.  Reasoning  along  the  lines  which  led  to 
finding  the  optimum  current  density  distribution  of  a  layer  coil  (see  Appendix  I)  leads  to 
the  following  optimization  condition  of  a  loop  coil  (compare  with  Eqs.  Ha  and  11b): 

^  _3  . 


+  3 

p  r^s 


=  constant  , 

(29a) 

1  sin  9^ 

p+  sin<?S, 

(29b) 

we  obtain: 

(31) 

Ja  P  ^  "  ' 

—  (loop  coil)  . 

(32) 

Substituting  this  value  in  Eqs.  27b  and  28b  we  obta 


\/2Tra, 


^opt 

Equation  29b  shows  that  the  current  density  distribution  does  not  depend  on  the  shape 
of  the  cross  sectional  area.  A,  of  the  coil  since  the  coordinates  of  the  contour  of  the  coil 
cross  section  do  not  appear  in  this  equation.  A  family  of  curves  can  be  drawn  to  repre¬ 
sent  the  i  =  constant  lines  which  yield  a  maximum  G  for  any  shape  of  coil  cross  section 
(see  Fig.  5a).  From  Eq.  31  it  follows  that  assumes  the  highest  values  if  the  coil 
cross  section  A  extends  from  points  as  close  as  possible  to  its  axis  of  symmetry.  This 
is  illustrated  by  the  following  examples  of  loop  coils  with  optimized  current  density 
distribution. 
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a.  Spherical  Coil  (Inside  Radius  Outside  Radius  a2,  see  Fig.  5b) 

A  spherical  coil  is  difficult  to  construct,  since  it  must  be  built  in  two  halves  in  order 
to  provide  access  to  the  inner  empty  sphere.  In  this  case: 

For  very  large  a  the  Fabry  factor  approaches  the  value  0.289. 

b.  Long  Cylindrical  Solenoid  (Inside  Radius  Outside  Radius 

c  ^/El  .  (34) 

20  y  a 

For  very  large  a  the  Fabry  factor  approaches  the  value  0.272.^^'^ 

c.  Cylindrical  Coil  with  Tapered  Ends  (Fig.  4) 

C  '  i„2V.-  4  Sin4*  ,  .  (35) 

°P*  10  V  4\2  /  2  16  y  a 


For  (/)=  7r/4: 


G  ,  =  0.262 

opt 


This  yields  =  0.213  for  a  =  3. 

As  shown  in  section  VIb  a  cylindrical  layer  coil  of  the  same  shape  optimized  in 
respect  to  current  density  distribution  has  a  Fabry  factor  =  0.204,  which  is  only 
slightly  smaller.  The  current  density  factor  J however,  increases  noticeably  from 
0.346  for  the  layer  coil  to  0.470  for  the  optimized  loop  coil.  This  shows  clearly  the 
merely  academic  interest  of  this  case. 

Finally  the  following  theorem  might  be  derived.  Consider  a  set  of  n  loop  coils  of  any 
cross  sections  (Fig.  6).  The  optimized  Fabry  factor  of  the  coil  k  is  (Eg.  31): 


opt  ""  jo 


y  ( sin  <!> 


The  optimized  Fabry  factor  of  the  whole  loop  arrangement  is; 


\j2TTay 
^opt  10 


( sin 


The  integration  has  to  be  performed  over  the  total  cross  sectional  area  A,  i.e.,  over  al 
the  incremental  areas  Expressed  in  somewhat  different  notation; 


^opt  ""  1^ 
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Fig.  6. 


Optimum  Current  Density  Distribution  of  a  Loop  Coil 
( />+  and  X"’"  Known  ). 

Figures  5a,  5b,  and  6. 

Therefore, 


(38) 


The  Fabry  factor  is  increased  by  the  addition  of  more  coils  to  the  set.  Its  highest 
value  is  achieved  when  the  whole  available  space  is  filled  with  windings.  This  is  in 
contrast  to  layer  coils.  For  these  the  Fabry  factor  decreases  when  a  certain  optimum 
winding  space  is  exceeded. 


CONCLUSION 

1.  General  expressions  for  the  Fabry  factor  of  a  “layer  coil*'  (for  which  coil  varia¬ 
bles  depend  on  the  distance  r  from  the  axis  only),  and  for  a  “loop  coil**  (coil  variables 
depend  on  r  and  on  axial  coordinate  z)  have  been  derived  (Eq.  7b  and  27b).  These  gen¬ 
eral  forms  are  convenient  for  the  computation  of  the  Fabry  factor  of  a  coil  of  any  shape 
and  any  distribution  of  current  density,  resistivity,  and  space  factor. 
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2.  In  order  to  produce  a  desired  field  strength  Hq  at  any  point  on  the  axis,  with  the 
minimum  power  P,  the  following  kinds  of  optimization  are  possible: 

a.  The  coil  shape  and  the  distribution  of  p  and  X  are  known;  the  problem  is  to  find 
the  optimum  current  density  distribution.  For  symmetrical  loop  coils  and  layer 
coils  with  variable  p  and  X  general  solutions  have  been  found  (Eq.  11b  and  29b). 
The  following  result  is  of  particular  interest:  A  desired  magnetic  field  strength 
in  the  center  of  a  symmetrical  layer  coil  is  achieved  with  the  minimum  power 
when  the  product  of  resistivity  p,  radial  distance  r  of  a  layer  from  the  axis,  dis¬ 
tance  of  the  layer  end  point  from  the  coil  center  s,  and  current  density  i  is 
constant. 

b.  The  current  density  distribution  of  a  layer  coil  Is  any  known  function  i  =  i(r), 
likewise  the  distributions  of  p  and  X  are  known;  find  the  shape  of  the  optimum 
coil.  A  general  solution  has  been  found  (Eq.  21).  Of  particular  interest  is  the 
current  density  distribution: 


In  this  case  the  optimum  coil  is  cylindrical  with  tapered  ends, 

c.  Simultaneously  optimizing  current  density  distribution  and  coil  shape  likewise 
leads  to  the  cylindrical  coil  with  tapered  ends  and  current  density  inversely  pro¬ 
portional  to  the  square  of  the  distance  from  the  axis. 

3.  For  any  magnet  coll,  the  current  density,  in  the  innermost  layer  (or  in  the  loop 
with  smallest  radius)  can  be  represented  by  the  product  of  a  dimensionless  factor,  the 
‘‘current  density  factor**  /,  and  a  scaling  factor  as  follows: 


'VpiV? 


(Ho) 


The  analogy  to  the  Fabry  factor  is  obvious.  General  expressions  for  the  current  density 
factors  of  layer  and  loop  coils  have  been  derived  (Eq.  14b  and  28b).  For  a  coil  of  any 
shape  with  constant  resistivity,  space  factor^  and  current  density,  the  current  density  fac¬ 
tor  J  is  the  reciprocal  value  of  the  square  root  of  the  volume  factor  (Eq.  15). 

4.  For  coils  with  optimum  current  density  distribution,  simple  general  expressions 
for  the  Fabry  factor  G  ,  have  been  derived  (Eq.  12  and  21).  The  current  density  factor 

'  opt 

of  a  layer  coil  becomes  in  this  case: 

j  =  — ! - (layer  coil)  ,  (17) 

and  that  of  a  loop  coil; 


/ 

Opt 
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5.  The  Fabry  factor,  optimized  with  respect  to  current  density,  of  n  loop  coils,  each 
having  the  optimized  Fabry  factors  is: 


opt 


opt 


(38) 


Numerical  examples  illustrate  the  basic  concepts  and  theorems  mentioned  above. 


APPENDIX  I 

Bitter^  showed  how  to  find  the  optimum  current  density  distribution  and  the  corre¬ 
sponding  values  of  the  Fabry  factor  of  colls  of  known  shape.  He  did  this  by  the  method 
of  indeterminate  multipliers,  and  the  same  mathematical  method  shall  be  applied  here  to 
the  optimization  of  the  current  density  z(r)  of  a  layer  coil  of  any  symmetrical  shape  and 
with  any  p{r)  and  X(r)  distribution. 

From  Eqs,  4  and  5, 


Att  r^2  z  . 

F  [i  (r)]  =  Hq  +  /xP  =  —  I  \i  —  dr  +  (I  Arr  I  pXi^rzdr  , 

^  ®  ^  ^  *^r=a  ^ 

p  being  an  arbitrary  constant. 


8F  [z(r)]  =  0 


(a) 

(b) 


is  satisfied  if: 

477"  2: 

- X — +p87TpXirz  =  0  .  (c) 

10  s 

Therefore: 

prsi  - - =  constant  =  iZ,  .  (d) 

20fi  ^1111 

This  is  Eq.  11a.  Note  thatX(r)  does  not  occur  in  this  optimization  condition  for  z(r). 

The  space  factor  A(r),  of  course,  influences  the  Fabry  factor,  as  shown  in  Eq,  12.  A 
similar  calculation  can  be  used  to  prove  Eq.  21. 


APPENDIX  II 

Equations  for  Unsymmetrical  Layer  Coils 

Figure  7  shows  a  layer  coil  with  the  z-axis  as  the  axis  of  rotational  symmetry.  The 
coil,  however,  is  not  symmetrical  to  any  plane  perpendicular  to  the  z-axis  (unsymmetrical 
layer  coll).  Referring  to  Fig.  7,  the  abbreviation, 

cos  (;6j|  -  cos  (?Sj  =  i/r  ,  (e) 

is  used.  The  magnetic  field  strength  at  the  point  0  assumes  the  form 


H 


0 


rr 
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Xil/i  dr  , 
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This  can  be  written  in  normalized  form: 


77  r^2  .  . 


The  ohmic  power  loss  becomes 


/: 


P  ~  2rr  I  p\r  (z  1 1  -  2 1)  dr 


=  27Tp^\jZ'2  r  -  z^){i*)^  dr  . 

From  these  equations  it  follows  that  the  Fabry  factor 

J*  ^  dr 


‘-k^R 


pVr(^,l 


and  the  current  density  factor 


/  = 


2^7  r ^  a 


.  .  -L.  O 
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The  current  density  distribution  for  which  G  is  a  maximum  can  be  found  by  the  method 
discussed  in  Appendix  I.  The  result  is 


p(^n  - 


i  -  constant  . 


(k) 


If  we  cons  ider  a  coil  cross  section  similar  to  that  shown  in  Fig.  1 ,  Eq.  (k)  applied  to  the 
point  T  -  a  ^  assumes  the  form  0/0,  which  evaluates  to 


p(g|i  - 

-A 


1  ' 

•  3  ' 

sm  cf)^ 


(I) 


or 


s  in 


01  ('^11  “  ^  I ) 


(m) 


This  value  substituted  in  Eq.  (i)  and  (j)  leads  to 


REFERENCES 

1.  A.  S.  Bishop,  Project  Sherwood.  The  U.S.  Program  in  Controlled  Fusion.  Addison- 
Wesley,  1958. 

2.  Proceedings  of  the  Symposium  on  Magnetic  Field  Design  in  Thermonuclear  Research^ 

held  in  Gatlinburg,  Tennessee,  December  11  and  12,  1958;  ORNL-2745. 

3.  The  ORNL  Thermonuclear  Program,  AEG  Research  and  Development  Report,  ORNL- 
2457. 

4a.  Ch.  Fabry,  L*Eclairage  Electrique,  XVII  No.  43,  133—41  (1898). 

4b.  Ch.  Fabry,  Journal  de  Physique,  IX  129-34  (1910), 

5.  J.  D,  Cockcroft,  Proc.  Roy.  Soc. ,  CCXXVIl  A654^  317-43  (1928). 

6.  F.  Bitter,  RSI,  VII  (1936),  479-87;  VIII  318-19  (1937);  IX  373-81  (1939). 

7.  Fabry,  ref  4a,  Figure  2. 

8.  W.  F.  Gauster,  Controlled  Thermonuclear  Conference  held  at  Washington,  D.C., 
February  3-5,  195S,  TID-7558,  pp  518-27. 

9.  Fabry,  ref  4a,  p  138. 

10.  Compare  for  instance  S.  S.  Attwood;  Electric  and  Magnetic  Fields,  John  Wiley  and 
Sons,  New  York. 

11.  Compare  Bitterns  computation  of  G  for  several  special  cases  (see  ref  6). 

12.  J.  C.  Maxwell,  A  Treatise  on  Electricity  and  Magnetism,  3d  Ed.  1892,  pp  360-4, 

13.  One  special  optimum  problem  treated  in  Fabry's  first  paper  4a  contains  an  error.  In 
Fabry's  second  paper  4b,  however,  all  optimum  values  indicated  are  correct. 


17 


